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Dephasing of Josephson qubits close to optimal points 
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Decoherence of Josephson qubits can be substantially reduced by tuning their parameters to optimal op- 
eration points, with only quadratic coupling to fluctuations. We analyze dephasing due to 1/f noise for a 
two-level system, detuned from an optimal point, i.e., the crossover to the linear-coupling regime, both for 
free induction decay and for spin-echo experiments. Influence of several noise sources is also discussed. 
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Superconducting nanocircuits are promising candi- 
dates for implementation of quantum-coherent two-level 
systems, building blocks of prospective quantum in- 
formation processing devices El- At present, co- 
herence of these circuits is limited mostly by low- 
frequency noise, often with 1/f power spectrum, such as 
background-charge fluctuations, variations of magnetic 
flux in superconducting loops, and of the critical current 
'of Josephson junctions pilEIE]- 

An efficient method of improving the coherence 
properties of qubits, the "optimal point strategy" , was 
suggested in the bias of a charge-phase qubit was 
tuned to an operation point where the linear coupling 
to the noise sources vanishes. As a result, the coherence 
time was extended by 2-3 orders of magnitude com- 
pared to earlier experiments. This approach was later 
applied to reduce the influence of bias-current fluctua- 
tions on a flux qubit 0- This method may be used in 
combination with the spin-echo- type techniques [SIEIIZIj 
inherited from NMR. 

The long coherence time in these devices allows to 
generate long-lived quantum-coherent oscillations and 
to study their decay laws in detail. Thus solid-state 
qubits may be used as a unique tool to gather infor- 
mation about the properties of the noise [HI EJ • Such 
studies El QUI allow one to investigate the origin and 
the microscopic mechanism of the noise and to design 
devices with even better coherence, necessary for large- 
scale quantum circuits. 

Analysis of decoherence at and close to an optimal 
point is important for superconducting quantum bits. 
Comparison of dephasing at an optimal point and away 
provides additional information on the nature and sta- 
tistical properties of the fluctuations (cf. Ref. Fur- 
thermore, it is relevan t since a two-level system may 
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be tuned away from degeneracy for quantum manipula- 
tions or readout. In other words, here we analyze the 
crossover between the optimal bias conditions and the 
linear-coupling regime far away from degeneracy. 

Recent experiments suggest that the low-frequency 
noise in Josephson circuits is produced by collections 
of bistable fluctuating systems, a well-known model of 
the flicker noise [TT]. For large collections with a suffi- 
ciently regular distribution of parameters, one expects 
the noise to obey Gaussian statistics due to the central 
limit theorem. For small number of fluctuators or sin- 
gular distributions, non-Gaussian effects may strongly 
influence decoherence El • 

Dephasing of qubits by long-correlated Gaussian 
noise at optimal points was studied earlier for 1/f noise 
spectrum in the cases of free-induction decay and 
echo m3 experiments (cf. also Ref. JH1)- Here we an- 
alyze decoherence in a system subject to Gaussian 1/f 
noise, from one or several independent sources, in the 
vicinity of an optimal point, such that both quadratic 
and linear coupling terms are relevant. 

We express the dephasing law of a system near an 
optimal point in terms of the decay law precisely at the 
optimal point. We show that the ratio of these two de- 
phasing laws is dominated by the low-frequency noise. 
This allows us to find the relevant laws for free- induction 
and spin-echo decay. We also generalize the analysis to 
the case of several noise sources. 

Below we assume that the infrared cutoff frequency 
of 1// noise is the lowest frequency scale in the prob- 
lem, and in particular, is much lower than 1/t for rel- 
evant times. In our analysis we exploit this fact shift- 
ing, under proper conditions, the spectral weights of the 
fluctuations between zero frequency and u)- lT . 
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We consider a two-level system (spin-1/2) with the 
Hamiltonian 

H = -^[e + Ae(X(t))]a z +H hath , (1) 

which describes coupling to a noise source via a fluctuat- 
ing quantity X(t). In turn, its dynamics is governed by 
the Hamiltonian of the noise source, Wbath- Physically, 
X(t) may represent, e.g., gate-voltage or magnetic-flux 
fluctuations in a Josephson qubit. We assume the 1// 
spectrum of its fluctuations in the relevant frequency 
range. 

Eq. (JTJ presents the Hamiltonian in the eigenbasis 
of the non-perturbed (Ae = 0) part. Only the fluctua- 
tions of the level splitting, i.e., the diagonal terms, are 
accounted for. Note that as far as the influence of the 
low-frequency noise (u <C £o) on the qubit is concerned, 
the general case reduces to Eq. in the adiabatic ap- 
proximation |17j . 

To characterize decoherence, we analyze the decay 
of the off-diagonal entry of the qubit 's density matrix 
in the eigenbasis. Averaging over the noise realizations, 
we find 

(a_(t))=Tr(a_p(t)) = (a_(0)>e ite op(t), 
Pit) = (T e*Jo 9(*) A£ ( x (*)) d *V 

where g(t) = 1 for the free induction decay. More gen- 
erally, the function g(t) accounts for modulations of the 
qubit-noise coupling. For example, in echo-type exper- 
iments \g(t)\ — 1 and the sign of g(t) is reversed every 
time when a 7r-pulse is applied. 

In the vicinity of an optimal point, keeping the 
leading-order terms in the expansion of Ae in X, we 
find 



A e (X(t)) = A 



Xit) + — 



(3) 



where we have combined the linear and quadratic terms 
to form a full square, by transferring a constant from 
£o for convenience. Here the constant D characterizes 
the offset from the optimal point and vanishes at the 
degeneracy. 

Dependence on the shift from the optimal point. One 
can represent the averaging in Eq. via a Gaus- 
sian functional integral ^Hl over X(t), of the type 

J e -Y T AY + D(eY+Y T t)+cD* dY _ Thig & relation 

between the decay law of coherence in the vicinity of 
the optimal point, P(t), and the dephasing law Po(t) at 
the optimal point (i.e., for D = 0): 



c 

FIGQ Diagrams that contribute to the dephasing at 
an optimal point (a) and close to it (b) 

with a D-independent function f(t). Naturally, this 
function depends on git) and on the noise power Sx{ui). 
Eq. (0J also follows from the diagrammatic analy- 
sis circular diagrams, considered in Ref. [Tl] . 
contribute to Po(t), and (the logarithm of) the sec- 
ond term is a combination of linear diagrams, with two 
linear-coupling vertices at the ends (Fig. 1), thus lead- 
ing to the D 2 dependence. Below we find f(t) in various 
situations. Note that in the far-detuned limit D — > oo 
only the short-time behavior of fit) is relevant for the 
dephasing law. 

Reduction to the dephasing law at the optimal point. 
Let us formally average the last equation over a 
Gaussian-distributed D with dispersion a 2 . Clearly, the 
average of the dephasing law Q over this distribution 
coincides with the dephasing induced by fluctuations 
with the modified noise power, S(ui) — S(ui) + ^a 2 S(ui). 

Physically, at the time scale t the low-frequency 
part of Xit) (i.e., the contribution of the frequencies 
u) <C 1/t) behaves as a static (u> = 0) random Gaus- 
sian quantity, and the dephasing law depends only on 
its dispersion, i.e., the low-frequency weight J du>S(u>). 
Technically, this follows, for instance, from the diagram- 
matic analysis, similar to that in Ref. |14j . In particu- 
lar, the addition of a new component at zero frequency 
is equivalent to the renormalization u>i T — > Lu' iT : 

A^ + -|ln(^ ir )| = -|ln(^ r )|, (5) 

where T = XX 2 , and Xf characterizes the strength of 
the 1/ f noise, 



Thus after averaging over D, one obtains 



P {t,u ir e-^ a ) 



[1 + 2\a 2 f(t)) 



1/2 



(6) 



(7) 



Comparing the subleading terms in the expansion in a 2 
at o — > 0, we find 2 > 



-m- 



* dlnPbif) (8) 

§rJjtW«m with a sharp low- 



P{t) = P (t)e- XD2 f^ 



(4) 



2 'For an arbitrary (non- 
frequency cutoff at a certain frequency u)\ T one finds f(t) 

7T d In Pp 

4AS(^ ir ) duj iT ' 
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Taking into account the dimensions of time t, noise 
power r, infrared w; r and ultraviolet lu c cutoff frequen- 
cies, one can rewrite Eq. (J7J) as 

_ £_ ( dhxPbjt) _ diiiP Q (t) _ dinP Q (t) \ 

I[ > 4T \ dint dlnT dlnu c J " 1 ' 

Eqs. (@J and JHJ allow one to find the dephasing close 
to an optimal point in terms of the dephasing law pre- 
cisely at this optimal point. 

Free induction decay, high ultraviolet cutoff. Short- 
time asymptotics of Po(t) were found earlier ^] I15| . 
For example, in the case of free induction decay and 
high ultraviolet cutoff, Lu c t 3> 1, one immediately finds 
from Eqs. |@J|, © and the results of Ref. 14 that 

-1/2 



P(t) =11- 



-iTt In — 

7T LU^-t 



/(*) 



it 



4 MTtln- 



(10) 



(11) 



These results are based on the short-time behavior of 
P (t) at Tt < 1 (and | ln(w ir t) > 1) [H]. They are 
dominated by the contribution of the low- frequency fluc- 
tuations (\uj\ <C 1/t). We show below that this contri- 
bution dominates at longer times as well; thus Eq. (jl 1|> 
(but not Eq. (fT~0|t ^ applies at longer times too. 

Indeed, we have seen that f(t) can be found from 
the sensitivity of Po(t) to w; r . Clearly, the contribution 
of high frequencies lo > 1/t is insensitive to the infrared 
cutoff. Before providing an accurate evaluation of f(t) 
let us remark that its low-frequency contribution can be 

estimated by treating the low-frequency part as static 

x 2 

noise with the same dispersion, er 2 = — ln^-r. Then, 
averaging the phase factor exp(i\(X + D/2) 2 ) over the 
Gaussian distribution oc cxp(— X 2 /2a 2 ), we find imme- 
diately the results (Ell) and {HJ ■ 

To calculate f(t) one has to evaluate the diagrams in 
Fig. lb. In contrast to the circular diagrams in Fig. la, 
the high-frequency contributions in Fig. 2b are sup- 
pressed since the incoming frequencies at the ends are 
zero, and at each vertex the typical frequency change 
is of order 1/t. Accurate evaluation demonstrates that 
this indeed suppresses all contributions apart from (|llf) . 
dominated by very low frequencies w«l/(. 

We remark that in the far-detuned limit D 2 — > oo 
the total decay Po(t) exp[— XD 2 f(t)] is dominated by the 
second factor and by the short-time asymptotics of f(t). 
In this limit one recovers, as expected, the decay laws 
for a linearly coupled reservoir. 

Free induction decay, low ultraviolet cutoff. For 
oj c t < 1, X(t) can be treated as static, and one eas- 
ily gets 



P (t)= I 1 - -itTln — 



-1/2 



/(*) 



it/ 4 



MtTm^- 1 



(12) 
(13) 



Spin echo decay near optimal points. For the case of 
free induction decay, the short-time asymptotics of Po{t) 
gave via Eq. (JHJ an expression for f(t) applicable at all 
relevant times. Similarly, one can find f(t) for a system 
subject to spin-echo pulses, using Eq. J3J) and the de- 
phasing laws at the optimal point in the corresponding 
cases (for one or several 7r-pulses and depending on 
the value and shape of the ultraviolet cutoff) . 

In particular, for N — 1 spin-echo pulses (N > 2) 
and a sufficiently high ultraviolet cutoff, uj c ^> N/t, one 
finds 



/(*) = 



N 



sr 



i + ^ (art)- in. 



(14) 



where Cjv is a constant of order 1 |15) . 

However, the derivation of Eq. ||SJ) for spin-echo de- 
cay requires stricter conditions: it implies that the de- 
phasing doesn't change significantly if the offset D from 
the optimal point is substituted by a slow oscillating 
offset with frequency uj^ -C 1/t. As a result of such 
a substitution the term A J g(t)D 2 dt = is substi- 
tuted by a non-zero quantity. This does not change 
the result as long as this quantity is much less than 
XD 2 f(t). An estimate of the inaccuracy of Eq. © leads 
to the constraint (twj r )t <C Nf(t) on its applicability 
to an echo-type experiment considered, which for rel- 
evant experimental parameters (lo- u . <C T) reduces to 
(ta;i r )t -€i A^/|rin(tw; r )|. For time scales of interest in 
spin-echo experiments this constraint is satisfied. 

Several fluctuating parameters. Typically, the be- 
havior of Josephson qubits is controlled by several pa- 
rameters, exhibiting low-frequency noise. The results 
derived above can be generalized to the situation with 
many sources of 1// noise in the vicinity of their optimal 
points. 

In the simplest case, various noise sources contribute 
independently: 



(15) 



Then the dephasing is a product of the partial contribu- 
tions, P(t) = Yli Pi(t), regardless of the noise spectra. 

However, in general the quadratic part of the per- 
turbation near the optimal point contains cross terms, 
XiXj. For instance, this is the case for the flux qubit 
in Ref. ^21 and two fluxes as control parameters. If all 
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noise sources have similar spectra (such that Si(ui) — 
const • Sj(u>)), one can reduce the problem to that for 
independent sources by simultaneous diagonalization of 
the quadratic perturbation and the correlation matrix 
(Xi(t)Xj{t')) (cf. Ref. PS])- For instance, this remark 
applies to several sources of 1// noise with the same in- 
frared cutoff (and, e.g., sufficiently high ultraviolet cut- 
offs, ui l c t > 1, which do not influence the dephasing law 
strongly) . 

In general, however, different sources of 1/f noise 
may have different infrared cutoffs. Still, this case can 
be reduced to the problem of a single noise source. In- 
deed, to treat this case, one can employ the approach 
used above in the derivation of JSJ): one adjusts all the 
infrared cutoff frequencies, w? r , to the same value (still, 
<C l/t) by shifting the rest of the low- frequency spectral 
weight to zero frequency. As a result, one arrives at the 
situation with several noise sources with the same spec- 
tra (up to an overall factor) and a static random field. 
At this point, one can average first over the dynamical 
fluctuations by re-diagonalization as above, and then 
over the static noise. Notice, that fluctuations with a 
low ultraviolet cutoff, <C l/t, may be treated as low- 
frequency fluctuations in the same spirit. 

In particular, the dephasing for a linear combination 
X(t) =J2 diXi(t) of 1// noise sources with high ultra- 
violet cutoffs can be obtained following the procedure 
described. One finds that the dephasing in this case co- 
incides with the dephasing due to a single source of 1/f 
noise with the power T = Yl |o;.i | 2 r^ and the infrared 
cutoff 

Lu iT = it (<4) Q * r</r . (16) 

In conclusion, we considered the dephasing of a two- 
level system close to an optimal operation point, in the 
crossover regime between the quadratic and the linear 
coupling to the noise source. We found the dependence 
on the shift from the optimal point, by relating the 
dephasing law to that at the optimal point with only 
quadratic coupling. In the case of 1// noise, we present 
an explicit expression for the ratio of two dephasing 
laws. Further, we analyzed the influence of several noise 
sources and demonstrated, how the problem can be re- 
duced to the case of a single source. 

We thank A. Shnirman for useful discussions. This 
work was supported by the Dynasty foundation and the 
grant MD-2177.2005.2. 
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